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1. Introduction 

An accurate determination of the thermodynamic properties of matter under extreme 
conditions of temperature and density is required for a correct description of the 
mechanical and thermal properties of many dense astrophysical bodies, including giant 
planets, low- mass stars (i.e., stars smaller than the Sun), and so-called compact stars 
(white and brown dwarfs and neutron stars). These objects are composed dominantly 
of ion-electron plasmas, where ions are strongly correlated and electrons are strongly or 
partially degenerate: classical Coulomb coupling parameter Fj = {ZicY /ksTai is large 
and electron density parameter = ai/ (aoZl^^) is less than unity (here ao = ^^/(mgC^) 
denotes the electronic Bohr radius, a-i = (3/47rni)^/^ the mean inter- ionic distance, 
the ion charge number, and ni the ion number density). The correct description of 
the structure and cooling of these astrophysical bodies thus requires the knowledge of 
the equation of state (EOS) and the transport properties of such dense plasmas. In 
this short review, we focus on the two extremes of this range of astrophysical objects 
in term of matter density: Jovian planets and neutron stars. As will be shown in the 
next sections, modern experiments and observations provide stringent constraints on the 
thermodynamic properties of dense matter under the physical conditions characteristic 
of these objects. 
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2. The equation of state of hydrogen and the structure of Jovian planets 
2.1. Hydrogen pressure dissociation and ionization 

Jupiter and Saturn are composed of about 70%-97% by mass of hydrogen and heUum. 
Temperatures and pressures range from T = 165 K and T = 135 K at P = 1 bar, 
respectively, at the surface, to T > 8000 K, P > 10 Mbar at the center. At pressures 
around P ~ 1-3 Mbar, corresponding to about 80% and 60% of the planet's radius, 
as measured from the planet center, for Jupiter and Saturn, respectively, hydrogen 
undergoes a transition from an insulating molecular phase to a conducting ionized 
plasma. The description of this transition, described as the pressure-ionization or 
metallization of hydrogen, has remained a challenging problem since the pioneering 
work of Wigner & Huntington ^ . Much experimental work has been devoted to this 
problem, but no conclusive result has been reached yet. Several high-pressure shock wave 
experiments have been conducted in order to probe the EOS of deuterium, the isotope of 
hydrogen, in the regime of pressure ionization. Gas gun shock compression experiments 
were generally limited to pressures below 1 Mbar j2], probing only the domain of 
molecular hydrogen. New techniques include laser-driven shock-wave experiments 
13 E], pulse-power compression experiments and convergent spherical shock wave 
experiments [3 IE] and can achieve pressures up to 5 Mbar in fluid deuterium at 
high temperature, exploring for the flrst time the regime of pressure-dissociation and 
ionization. These recent experiments give different results at P > 1 Mbar, however, and 
this controversy needs to be settled before a robust comparison between experiment and 
theory can be made in the very domain of hydrogen pressure ionization. 

On the theoretical front, a lot of effort has been devoted to describing the pressure 
ionization of hydrogen. The EOS commonly used for modeling Jovian planet interiors 
is the Saumon-Chabrier-Van Horn (SCVH) EOS jHl GDI HI] wich includes a detailed 
description of the partial ionization regime. This EOS reproduces the Hugoniot data 
of Nellis et al [21 but yields temperatures about 30% higher than the gas reshock data, 
indicating insufficient D2 dissociation 12j. A slightly revised version jTSj recovers the 
gas gun reshock temperature data as well as the laser-driven shock wave results [3], 
with a maximum compression factor of p/po — 6, where po = 0.17 g cm~^ is the 
initial density of liquid deuterium at 20 K. On the other hand, the earlier SESAME 
EOS Jl], based on a similar formalism, predicts a smaller compression factor, with 
p/po ~ 4, in general agreement with all the other recent shock wave experiments. Ah 
initio approaches for the description of dense hydrogen include path integral Monte 
Carlo (PIMC) ITT)| IT7j and quantum Molecular Dynamics (QMD) simulations. The 
latter combine molecular dynamics (MD) and Density Functional Theory (DFT) to take 
into account the quantum nature of the electrons [T8 t ll9 | l2( H l21j. The relevance of earlier 
MD-DFT calculations was questioned on the basis that these simulations were unable to 
reproduce data from gas-gun experiments ^Hl- This problem has been solved with more 
accurate simulations [ini EUl EI] • Although an ah initio approach is more satisfactory 
than the phenomenological approach based on effective potentials, in practice these 
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simulations also rely on approximations, such as the handling of the so-called sign 
problem for the antisymmetrization of the fermion wave functions, or the calculation of 
the electron functional density itself (in particular the exchange and correlation effects), 
or the use of effective pseudo-potentials of restricted validity, not mentioning finite size 
effects. Moreover, these simulations are too computationally intensive for the calculation 
of an EOS covering several orders of magnitude in density and temperature, as necessary 
for the description of the structure and evolution of astrophysical bodies. 

Figure Q compares experimental and theoretical Hugoniots in the P-p and P-T 
planes. The disagreement between the laser-driven experiments and the other techniques 
is illustrated in the P-p diagram. Whereas the SCVH EOS achieves a maximum 
compression similar to the laser-driven data, all the other models predict compression 
factors in the P-p plane in agreement with the more recent data. The MD-DFT 
results, however, predict temperatures for the second shock significantly larger than 
the experimental results jT^j. Even though the experimental double-shock temperature 
may be underestimated due to unquantified thermal conduction into the window upon 
shock reflection, and thus represents a lower limit on the reshock temperatures, the 
disagreement in the T-V plane is significant. As noted previously, the degree of 
molecular dissociation has a significant influence on the thermodynamic properties of 
the fluid and insufficient dissociation in the simulations may result in overestimates of 
the temperature. It has been suggested that the LDA/GGA approximations used in 
MD-DFT underestimate the dissociation energy of D2 This would lead to even 

less dissociation. The fact that compression along the experimental Hugoniot remains 
small thus suggests compensating effects in the case of hydrogen. More recent, improved 
simulations |;2lj, however, seem to partly solve this discrepancy and to produce reshock 
temperatures in better agreement with the experimental results. Peak compression 
in the modern MD-DFT simulations occurs in the ~ 0.2-0.5 Mbar range around a 
dissociation fraction of ~ 50%. 

The differences in the behaviour of hydrogen at high density and temperature 
illustrated by the various results displayed in Fig. bear important consequences for 
the structure and evolution of our Jovian planets. These differences must be correctly 
understood before the description of hydrogen pressure dissociation and ionization 
stands on firm grounds. As noted by Boriskov et al |2ni, all the recent experiments 
agree quite well in terms of the shock speed Us versus the particle velocity Wp, almost 
within their respective error bars. Error bars and differences in {us,Up) are amplified 
in a P-p diagram by a factor of (p/po — !)• These are challenging experiments as the 
differences seen in panel 1 of Fig. 1 arise from differences in Ug and Up of less than 3%. 
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Figure 1. Experimental sliock (P, p, T) data and tlieoretical Hugoniots of deuterium. 
Sources of data are: Gas gun et EHE], Z machine NOVA 011] and CSSW 
[3 Eg. Curves show Hugoniots computed from the EOSs of SCVH [TTl, SESAME 
mi, PIMC HSI, and MD-DFT [201 ■ 



2.2. Jupiter and Saturn interiors 

The rapid rotation of Jovian planets induces a nonspherical gravitational field that can 
be expanded in Legendre polynomials PnicosO): 



GM 



n 



yM = -^[^-Y.[-f) JnPn{cO,e)[ (1) 

n=l ^ ' 

where M and Req denote respectively the planet mass and equatorial radius, and the 
Jn are the gravitational moments: 

Jn = -j^ 1^ r-P„(cos e)p{r', 6) rfV. (2) 

Because of north-south symmetry, the moments of odd order are null. The first 
three nonvanishing moments, J2, J4 and Jg have been measured with high accuracy for 
both planets during spacecraft flyby missions. Combined with the planet mass, radius 
and rotation period, these provide integral constraints on the density profile of the 
planet, p{r,9), to be compared with the corresponding values from a structure model 
obtained for a self-gravitating and rotating fluid body in hydrostatic equilibrium. The 
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EOS provides the -P(p) relation needed to close the system of equations. The structure 
of the H/He envelopes of giant planets is fixed by the specific entropy determined from 
observations at their surface. The very high efficiency of convection in the interior of 
these objects leads to nearly adiabatic interior profiles. The structure of the planet is 
thus determined by the choice of the hydrogen and to a lesser extent by the helium 
EOS used in the models. A detailed study of the infiuence of the EOS of hydrogen on 
the structure and evolution of Jupiter and Saturn has been conducted recently |24j . 
Fortunately, some shock wave experiments overlap Jupiter's and Saturn's adiabats. 
Figure El displays Jupiter (J) and Saturn (S) adiabats for hydrogen calculated with 
the SCVH EOS and the first and second shock Hugoniots calculated with the SESAME 
EOS and illustrates relative differences in density between Jupiter adiabats computed 
with these two EOSs. As demonstrated by Saumon & Guillot j2lj, the small (< 5%) 
difference on the (P, p) relation along the adiabat between the two EOSs, representative 
of the two sets of experimental results, is large enough to affect appreciably the interior 
structure of the models. Note that the SESAME D2 Hugoniot at low density is somewhat 
stiffer than the gas gun experiments [2j and does not recover the ideal H2 entropy at low 
temperature and density. No model of Jupiter could be obtained with this EOS 24j. A 
slightly modified SESAME EOS, which does recover the H2 entropy at low temperature 
and density, yields Jupiter models with a very small core mass, Mcoro ~ 1 ( M0 is 
the mass of the Earth) and a mass Mz ~ 33 M0 of heavy elements {Z > 2) mixed in 
the H/He envelope. The SCVH EOS yields models with M^^re ~0-6 M® and Mz ~ 15- 
26 M^. Models of Saturn are less sensitive to the EOS differences, since only ~ 70% of 
its mass lies at P > 1 Mbar, compared to 91% for Jupiter. Models computed with the 
SCVH and the modified SESAME EOS have Mcore = 10-21 M® and Mz = 20-27 M® 
and 16-29 M^, respectively. As seen in Fig. |2l the temperature along the adiabat is more 
sensitive to the choice of the EOS. This affects the thermal energy content of the planet 
and thus its cooling rate and evolution. Equations of state which are adjusted to fit the 
deuterium reshock temperature measurements [23] lead to models that take ~ 3 Gyr for 
Jupiter to cool to its present state. Even when considering uncertainties in the models, 
or considering the possibihty of a H/He phase separation, such a short cooling age is 
unlikely to be reconciled with the age of the solar system. This astrophysical constraint 
suggests that the reshock temperature data are too low. 

2.3. Helium equation of state. Plasma phase transition 

The planet interior models are also affected, to a lesser extent, by the uncertainties of the 
helium EOS. A model EOS for helium at high density, covering the regime of pressure 
ionization, has been developed recently by Winisdoerffer and Chabrier [221 • This 
EOS, based on effective interaction potentials between He, He+, He"'"+ and e~ species, 
reproduces adequately experimental Hugoniot and sound speed measurements up to ~ 1 
Mbar. In this model, pressure ionization is predicted to occur directly from He to He"''^. 
Because of the uncertainties in the treatment of the interactions at high density, however. 
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Figure 2. Adiabats for hydrogen in P-p and P-T planes. The curves labeled "S" and 
"J" show the SCVH-interpolated EOS adiabats of Saturn and Jupiter, determined 
by T = 136 K and T = 170 K at P = 1 bar, respectively. The first and second- 
shock Hugoniots calculated with the SESAME EOS are shown by the heavy solid line 
labeled " H" . The light solid curve (right-hand scale) shows the difference between 
Jupiter adiabats calculated with the SESAME EOS relative to the SCVH-interpolated 
EOS. 

the predicted ionization density ranges from a few to ~ 10 g cm~^. Comparison of the 
model predictions with available measurements of electrical conductivity of helium at 
high density j2Zl EH] is under way. 

The pressure ionization and metallization of hydrogen have been predicted to occur 
through a first order phase transition, the so-called plasma phase transition (PPT) 
[H ESI EOl ISD CHI I12]. Nearly all of these PPT calculations are based on chemical EOS 
models. Such models are based on a Helmholtz free energy that includes contributions 
from 1) neutral particles (atoms and molecules), 2) a fully ionized plasma, and 3) usually 
a coupling between the two. It is well-known that realistic fully ionized plasma models 
become thermodynamically unstable at low temperatures and moderate densities. This 
is analogous to the behavior of expanded metals at T = that display a region where 
dP/dp < and even P < [32] • This behavior of the fully ionized plasma model is 
formally a first order phase transition and reflects the formation of bound states in the 
real system. In other words, the chemical models have a first order phase transition 
built in from the onset, and this phase transition coincides, not surprisingly, with the 
regime of pressure ionization. This represents a common flaw in this type of models 
and it follows that their prediction of a PPT in hydrogen is not credible. This is further 
supported by a detailed study of two of these models [Hll 1121 • C)n the other hand, 
recent ah initio simulations find a sharp (6 ± 2%) volume discontinuity at constant 
pressure [HBl 121] or dPj dT < at constant volume EH ES] , a feature consistent with 
the existence of a first order phase transition. At the same time, the pair correlation 
function exhibits a drastic change from a molecular to an atomic state with a metallic 
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character (finite density of electronic states at tlie Fermi level). These transitions are 
found to occur in the ~ 0.5-1.25 Mbar and ~ 1500-3000 K temperature range. While 
these results are suggestive, a systematic exploration of this part of the phase diagram 
remains to be done. Note that a first order structural transition for H2 at T = is 
predicted to occur at a pressure P > 4.0 Mbar, from DFT calculations based on exact 
exchange calculations |22]- There is so far no published experimental evidence for the 
PPT but it cannot yet be ruled out. Given the difficulty of modeling this region of the 
phase diagram of hydrogen, only experiments can ultimately establish whether a PPT 
exists or not. 

3. Dense matter in strong magnetic fields. Neutron star structure and 
cooling 

Neutron stars (NS) consist of a nucleon core and envelopes composed of Coulomb 
plasmas. Cooling rates of these stars are determined by the heat capacity and neutrino 
emission processes in their cores and by heat transport in the envelopes. For the neutrino 
emission, most important are so called direct Urea (Durca) processes (beta-decay and 
beta-capture) and modified Urea (Murca) processes (the same but with participation 
of an additional nucleon, which helps to fulfill momentum conservation). The Murca 
processes are less powerful, but they work in every sufficiently hot NS. At contrast, 
the most powerful Durca processes operate only if the proton fraction in the core is 
large enough (otherwise the momentum conservation condition for degenerate nucleons 
cannot be satisfied). Some models of nuclear matter predict that a NS with relatively 
high mass should have a sufficient proton fraction at the stellar center for the Durca 
processes to occur. Such stars should cool faster, which opens a possibility to check the 
EOS of superdense matter through observations. The cooling rates are also strongly 
affected by nucleon superfiuidity (see ^| for review and references). 

Most NSs have magnetic fields B ~ 10^^-10^^ G, whereas some (so-called 
magnetars) are thought to have fields as high as ~ 10^^-10^^ G. The photosphere of a NS 
is characterized by temperatures Tg ~ 10^-10'' K (depending on the age t and mass M 
of the star) and densities p ^ 10~^-10^ g cm~^ (depending on T and B). Traditionally 
the NS crust is assumed to be composed of iron. However, the outer layers, including 
the atmosphere, can be composed of light elements (H, He, C) accreted on top of the 
iron layer. Therefore the determination of the temperature profiles and emitted spectra 
of NSs requires an accurate description of the formation of bound states and pressure 
ionization of these elements in a strong magnetic field. 

The quantum-mechanical properties of free charged particles and bound species 
(hydrogen atoms and molecules) are strongly modified by the magnetic field, which 
thereby affects the thermodynamic properties of the plasma [331 EHl- The transverse 
motion of electrons in a magnetic field is quantized into Landau levels. The energy of the 
nth Landau level of the electron (without the rest energy) is mgC^ (vT+26n — 1), which 
becomes hwcn in the non-relativistic limit, where hwc = hcB/nieC = 11.577-Bi2 keV, 
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is the electron cyclotron energy, b = hwc/meC^ = i?i2/44.14 is the field strength in 
the relativistic units, and B12 = -8/(10^^ G) is a typical magnetic-field scale for NS 
conditions. The atomic unit for the magnetic-field strength is set by huc = e^/ao, i.e., 
Bq = {rriec/he) x (e^/ao) = 2.35 x 10^ G. It is convenient to define a dimensionless 
magnetic-field strength 7 = B/Bq = b/ajr, where aj is the fine structure constant. 

For 7 ^ 1, as encountered in NSs, the ground-state atomic and molecular binding 
energies increase as ~ In^ 7. The H atom in a strong magnetic field is compressed in 
the transverse directions to the radius ~ Om, where 

= {hc/eBf/^ = 7-^/^ao = 2.56 x IQ-^^B^^^^^ cm (3) 

is the quantum magnetic length, which becomes the natural length unit. The increase 
of binding energies and decrease of sizes lead to a significant increase of the fraction of 
non-ionized atoms in the plasma at the photospheric densities (which are the higher the 
stronger the magnetic field). For example, at T = 10® K and B = 10^^ G, the typical 
density is p ~ 1 g cm~^, and there are > 1% of atoms in the H atmosphere. Because of 
the alignment of the electron spins antiparallel to the field, two atoms in their ground 
state (m = 0) do not bind together, because of the Pauli exclusion principle. One of 
the two H atoms has to be excited in the m = —1 state to form the ground state of the 
H2 molecule [SEj- Another important effect is that thermal motion of atoms across the 
field strongly modifies their binding energies and radiative transition rates. As shown 
in j37[ EHl inn] , the allowance for partial ionization and thermal motion is crucial for 
neutron-star atmosphere modelling. 

As long as T < huJkB = 1.343 x 10® Bu K and p < ~ 7.1 x 10^ B%'^ g cm'^ 
the electron cyclotron energy hwc exceeds both the thermal energy k-^T and the electron 
Fermi energy /cb^f, so that the field is strongly quantizing (e.g., In this case, typical 

for the NS photospheres, the electron spins are aligned antiparallel to the field. The 
electron Fermi energy decreases, therefore the onset of degeneracy is shifted to higher 
densities (slightly below p^). Proton motion is also quantized by the magnetic field, 
but the corresponding cyclotron energy is much smaller, huJcp = ^c'^e/''^p- 

A model which describes the thermodynamics of an interacting (H2, H, H+, e~) 
plasma in a strong magnetic field was constructed in |nZ|. On the base of this model, 
the EOS for magnetized H atmospheres of NSs, as well as their opacities, were explored 
and tabulated in ISHHSHI- 

Landau quantization of electron orbits affects not only the EOS and the radiative 
opacities, but also the heat conduction in the surface layers (see jHH] and references 
therein). The EOS of strongly magnetized, partially ionized hydrogen plasma as well as 
the electron conductivities and radiative opacities in neutron star magnetized envelopes 
were used in [101 to calculate the thermal structure and coohng of superfiuid NSs with 
accreted envelopes in the presence of strong dipole magnetic fields. In jlUj (see also jH] 
and references therein), the effect of neutron superfiuidity in the NS inner crust was also 
examined. The account of the effects of accreted matter, magnetic field and neutron 
superfluidity alters significantly the NS cooling. 
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Figure 3. Effective surface temperature (as seen by distant observer, T^) versus 
NS age t for assumed NS mass M = 1.3 and 1.5 solar masses. The dots with error 
bars show estimates of NS ages and effective temperatures from various observations; 
the dots with errors indicate observational upper limits. Left: Cooling of NSs with 
different relative masses AM/M of accreted (H-He-C) matter (values of log AAf/Af are 
indicated near the curves) Solid curves refer to non-accreted (Fe) iron envelope of the 
star. Right: Cooling of NSs with iron envelope for different magnetic field strengths 
(log_B in Gauss). 

Figure El displays theoretical cooling curves of NSs with (lower curves, M = 
1.5 Mq) or without (upper curves, M = 1.3 Mq) Durca processes in the core, with 
or without accreted envelopes, and with magnetic field of different strengths, compared 
to the estimates of the effective temperature obtained from observations (see jSj for 
references). As seen in the figure, the presence of a light-element (accreted) envelope 
increases Tg at the early cooling stage (t < 10^ yr), and as a result the thermal energy 
becomes exhausted sooner. The magnetar-like magnetic field B > 10^^ G acts in a 
similar way, whereas a weaker field almost does not affect the cooling. 

For simplicity, in Fig. IHlwe neglect the effects of superfluidity. Their discussion can 
be found in [SH EOl E] . 

4. Conclusions 

In this brief review, we considered the description of the thermodynamic properties of 
dense Coulomb matter in two specific astrophysical contexts, Jovian planets and neutron 
stars. The description of the pressure ionization of hydrogen and other elements, as well 
as the presence of strong magnetic fields, play an important role in determining the 
mechanical and thermal properties and the evolution of these objects. Models including 
these complex effects can successfully explain a variety of observations. On the other 
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hand, modem experiments and/or observations can enable us to discriminate between 
various EOS models in planet interiors and lead to a better determination of masses of 
accreted envelopes, surface magnetic fields and eventually the EOS of superdense matter 
in neutron stars. 
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